Quantum fluctuation theorem in an interacting setup: point contacts in fractional 

quantum Hall edge state devices 
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We verify the validity of the Cohen- Gallavotti fluctuation theorem for the strongly correlated 
problem of charge transfer through an impurity in a chiral Luttinger liquid, which is realizable 
experimentally as a quantum point contact in a fractional quantum Hall edge state device. This is 
accomplished via the development of an analytical method to calculate the full counting statistics 
(FCS) of the problem in all the parameter regimes involving the temperature, the Hall voltage, and 
the gate voltage. 
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Fluctuation theorems (FT), which capture universal 
properties of systems far from equilibrium, have been the 
subject of intense activity in the past fifteen years, both 
in the classical and in the quantum case One of 
the most celebrated examples is the Cohen-Gallavotti 
FT (CGFT) @, H about the ratio of the probability of 
having a time averaged entropy production s t in a steady, 
noncquilibrium state, take the value s and the probability 
that it takes the opposite value — s. The CGFT says that 
this ratio is a simple exponential: 



P(st = s) 
p(s t = -s) 



(1) 



Noncquilibrium fluctuations arc most naturally observed 
in open quantum systems exchanging energy or matter 
with reservoirs. Nanoscale electronic devices, where fluc- 
tuations can be resolved at the single electron level, pro- 
vide an ideal setup to investigate experimentally and the- 
oretically the CGFT. It is indeed easy to relate (JIJ) to a 
symmetry of the full counting statistics (FCS) generating 
function. Imagine for instance electronic transport be- 
tween two leads at the same (inverse) temperature /3 but 
different chemical potentials, coupled by a certain tun- 
neling hamiltonian. The generating function of charge 
transferred (for more precise definition, see below) is then 
expected to obey, in the large time limit [l[ 



x(k) = X (PV - k) 



(2) 



where V is the difference of potential, and we have set the 
electronic charge equal to unity. The verification that ((SJ 
holds in the noninteracting case for the so-called Lesovik- 
Levitov (LL) formula (3| is a corner stone of the rapidly 
developing field of electron counting statistics [f| . This 
formula reads 



lnx(/c) = t / de n(e) In 



-r(e) 



cosh(ft — 
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where t is the waiting time, n(e) is the density of states 
per unit length (or time, after setting the value of sound 
velocity to unity) at energy e, and r(e) is the modu- 
lus square of a scattering amplitude associated with the 
tunneling process. Expansion in powers of k gives rise 
to highly non trivial sum rules between the cumulants of 
order p and all cumulants of higher order. 

It certainly seems desirable to investigate relation fl2J) 
in the presence of interactions. Unfortunately, this re- 
quires, on the theoretical side, tackling interacting, out 
of equilibrium, strongly fluctuating open quantum sys- 
tems, which is in general a daunting task even for clas- 
sical systems @. We can nevertheless rely on a lot of 
recent progress. The idea of using the Bcthc ansatz to 
describe transport through quantum impurities has been 
pursued in several works |7|, [8|, leading, in some cases, to 
the calculation of full I — V characteristics Q , the second 
(noise) [Io[ and third cumulants of the current, and 
T = full counting statistics [l2[. While the approach 
led to some controversies, it has been confirmed repeat- 
edly by perturbative expansions 12 j_ 13(] and thorough 
numerical t-DMRG calculations [9. In fact, a very 
recent paper reports verification of the formulas in |12| 
for the FCS at T = [l5| over the whole crossover region 
between weak and strong effective couplings. 

Building on this progress, we report in this paper the 
calculation of FCS at non zero temperature in an inter- 
acting problem, together with the verification of the CG 
symmetry @. 

The basic system we consider is the ID chiral Luttinger 
liquid (CLL) with an impurity, which describes in partic- 
ular tunnelling between fractional quantum Hall (FQH) 
edges [3], or tunneling in the interacting resonant level 
model (IRLM) at the self-dual point 0] . 

The bulk Hamiltonian can be written in terms of 
right/left (R,L) moving charge current densities, 



9 



dx (jl 



fn) 
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where we fix h — 1 and where g is the Luttinger liq- 
uid interaction parameter [17j. In the case of a FQH 
device g = v when \ jv is an odd integer fl6j . The origi- 
nal fermions can be rewritten as exponents of the phase 
fields jlm, = TdxfpLM./^, so that the (local at x = 0) 
backscattering term (or intercdge tunnelling) is given by 



H B s = Ass cos [</>l(0) - <Pb{ 



(5) 



where Xbs is the respective amplitude. It has been shown 
in 0], that H + Hbs can be mapped (after a folding 
and a transformation into an even-odd basis) onto the 
boundary sine-Gordon (BSG) model, which is intcgrable 

The strategy to calculate properties in the steady state 
relics on a scattering description of the transport pro- 
cess. This is partly related to the Landauer-Biittiker 
approach, although the latter was developed for free the- 
ories, while we are dealing here with a (strongly) inter- 
acting situation. The point is that there exists a ba- 
sis of quasiparticles for the Hilbert space where: (a) the 
initial reservoir-like boundary conditions can be readily 
applied; and (6) scattering through the impurity is sim- 
ply described in a way that allows control of the charge 
flowing through the system. More precisely, this ba- 
sis, at 1/v odd integer - to which we restrict now - 
is made of kinks, antikinks (denoted by subscripts ±), 
and j = 1 . . .1/v — 2 breathers. All of them are mass- 
less and have the dispersion relations e±{9) = e e /2 and 
ej(6) = sin[7rj/2(l/^ — 1)] e e , where 9 is referred to as 
rapidity and parameterizes the particle momenta. While 
the breathers are neutral, the (anti)kinks carry charge (in 
units of the electron charge) Q = =pl, where by charge 
we mean the difference between the original charges in 
the L and R channels of (j4]). 

Of course, these quasiparticles are very complicated 
objects in terms of the bare electrons, but since all one 
wants to follow is the distribution of the charge, this does 
not matter. The point is that, because of integrability, 
the quasiparticles scatter without particle production, 
and the process in the large waiting time limit is entirely 
encoded in a scattering matrix, with all the one particle 
- diagonal as well as off-diagonal - elements given in [?| . 
We define the FCS generating function in the way dis- 
cussed for instance in [l9[ , by directly using the charge Q 
(instead of the currents). This charge enjoys both ther- 
mal and quantum fluctuations. We assume the system 
is initially prepared with decoupled channels at potential 
±V/2, and we are interested in the large time limit, where 
the system settles in a stationary process. A key feature 
of the integrable quasiparticle basis is that the DC com- 
ponent of the current acts diagonally on the multiparticle 
states, and the contributions to the FCS of different ra- 



pidities in the stationary limit simply factorize [2l| • The 
argument in [l9j can then be straightforwardly general- 
ized at every rapidity. Like in the free fermion case, one 



obtains in the end an expression for the FCS as the aver- 
age of a certain, k dependent expression calculated with 
respect to an equilibrium distribution in the reservoirs. 
Going through the mapping of the initial problem onto 
the BSG model one ends up with 

xM=(n{ i+T ( e )[( eK - i )^( i - ? ?+) 

e 

+ ( e - K -l> ?+ (l- ?/ _)]}), (6) 

where the r\± are occupation numbers (equal to zero or 
one) for single particle kink and antikink states of mo- 
mentum (energy) e in the equilibrium (and unfolded) in- 
tcgrable model. The brackets denote average over the 
equilibrium distribution determined by the potential V, 
and r(e) is the probability that a kink (antikink) goes 
through without being scattered into an antikink (kink). 

Hidden in the apparent simplicity of Eq. ^ is a funda- 
mental subtlety: the quasiparticles are not free, and we 
must think of them as a gas of particles whose wave func- 
tions satisfy a rapidity dependent statistics. This is the 
price to pay for the simple description of the tunneling 
physics in terms of single particle scattering. The prop- 
erties of bulk wave functions depend on the bulk sine- 
Gordon scattering matrix, whose elements in the massless 
limit are also given in ffy. The average in Eq. ([6]) there- 
fore can only be done using the thermodynamic Bethe 
ansatz technique. This means that the different levels e's 
in the product are not independent, but are solutions of 
a system of coupled Bethe ansatz equations. Like in the 
nonintcracting systems one can introduce the level densi- 
ties (per unit length and rapidity) rii(9) and filling frac- 
tions fi(9), i = ±, 1, . . . 1/v — 2, which can be combined 
to give the density of occupied states Pi(9) = rii(9)fi(9). 
But in contrast with the noninteracting case, the level 
density at rapidity 9 depends on the filling fractions at 
all other rapidities. 

Formula ([6]) can immediately be used to recover known 
expressions for the current and the DC noise. In the for- 
mer case for instance, expanding to first order in k gives 
readily / = I(V,T B ) = J d9 T {9){{P + - P-){9). Here 
(P±) denotes the average densities of kinks and antikinks 
in equilibrium on the line (without scattering) at temper- 
ature T and potential V. The scale Tg is a characteristic 
energy scale for the tunneling process, formally analog to 
the Kondo temperature in the Kondo problem. Its de- 
pendence on the microscopic tunneling amplitude Xbs m 
Eq. ([5]) is not universal; by scaling one has Tb oc X^g v . 
The physics of this expression is well known, and de- 
scribes the crossover between a regime where essentially 
no current is backscattered and / m vV to a regime where 
most current is backscattered. While the known expres- 
sion for the DC noise [lo| can be recovered as well, Eq. (J6j) 
gave readily access to the third cumulant, which was an- 
alyzed in 
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The calculation of the full FCS from Eq. ([6]) is a chal- 
lenging problem, and we will only give the main ideas 
here, postponing details to further publication. The first 
step is to write the generic term in the product in Eq. 
as an exponential using the simple identity 

1 + r [(e K - l)r?_(l - 77+) + (e- K - l)r?+(l - »/_)] 

= exp[m7_(l -??+)+ (1 -!?_)] (7) 

with a(f3) = In [1 + r(e ±K — 1)]. The next step is to go 
to a continuum limit where we replace discrete sets of al- 
lowed momenta by continuously varying rapidities, with 
associated densities. Note that in this replacement, the 
effective length of the system is in fact the waiting time t 
(where we have set vf = 1)- While a naive replacement 
?7± — > f± was used in it is not sufficient beyond the 
third cumulant. The point is that, in a continuum de- 
scription, with a certain number of allowed levels n± (9)d9 
for the ± species, and a certain fraction f±(9) of these 
levels occupied, the 'overlap' r/ + (8 i )r]_(9 i ) for the (still 
discrete) Qi € [0,9 + d9] can still strongly fluctuate. It is 
thus necessary to perform the sum over all the values of 
this overlap first, which can be done, in the large waiting 
time limit t — > oo, by a saddle point technique [22|. Once 
this is done, we can transform the average into a func- 
tional integral over the densities [20j , taking into account 
the usual TBA entropic contributions. The end result is 
of the form 



X = 



Z 1 J T>[n±,m; f±,fi] exp (s - ^ + S^j , 



(8) 



where Z is the expression in the numerator without 
S K , the energy has the form E = tj^if d9ei(9)Pi(9) 
with e±(9) — e e /2 + V/2 (the breather energies be- 
ing unaffected by the voltage) and the entropy S = 
-t J2i J d9 ni {9) [fi In ft + (1 - ft) ln(l - /<)] (9). ' The 
last term is of the form 



given by 



A(/+, /_) = (<* + /?)/-/+ + In 



(1 -/+)(!-/_) 



E U ln 



j=± 



1 - /+ - /- 



/) 



a -fm-f) 



,(io) 



where / satisfies the equation 



oa+[j = /(i _/ + -/_+/) 
(/+ - /)(/_ - /) 



111) 



and represents the distribution of states populated by 
both kinks and antikinks. 

Technically the evaluation of the functional integrals 
(JSJ) is done in several stages. The saddle point conditions 
with respect to filling fractions produce 1/v equations 
for the required 2/v variables n± j rii and P+,Pj. The 
Bethe ansatz equations (see e. g. [7J for further details) 
then yield the relations between densities of states and 
densities of occupied states Pj and are exactly the same 
as for the evaluation of conventional observables, 



E 

3 



(12) 



where i,j refer to (anti)kinks as well as breathers and 
&ij (9) are bulk scattering matrices given in Q . The den- 
sities n± are identical. 

We performed this program for the nontrivial case of 
filling fraction v = 1/3, where there is only one breather 
degree of freedom. The saddle point equations have then 
the following form 



Xi(6) = J d6'$i 



-6')Yi{9') + § 1+ {9-9')Y Q (6') 



S K = t[n(9)d9[af-(l-f + ) + (3f + (l-f-) + A(f + ,f_)] = j ' & - 9')Y l {9>) + $++(0 - 9')Y (9') 

(9) X-(9) = X+(9), (13) 
where n{9) = n±(9), and A takes into account the fluc- 
tuations of overlaps between ± species in Eq. (|6]). It is where Y± = ln(l — fx), Xi = —e\/T + ^(f^ 1 — 1), and 



X± 
Y 



= (3, a 



e± T V/2 , . (I -/+-/_ + / 
— \- In 



ln(l -/+-/_+/) + 



f±-f 

(!-/--/+)(/+ + /-) 
1 -/+-/-+/ 



J^^L + fJh - V (i-°L\ (i4) 

h - f df± (1 -/--/+ + f)(f± - /) V df ± > 




and represent a closed equation system with respect to filling fractions /j. They are conveniently solved for nu- 
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FIG. 1: (Color online) FCS generating function normalized 
to 2tTs/h for V/Tb = 2 and different temperatures. 



H 4 

S 

~ 2 







\ 










V/T = 

' B 

V/T -- 

' B 

V/T = 

' B 

V/T = 

' B 


=0.75 
=1.5 
=3 
-A 




i 




V 

\ 


\ 

\ 

\ \ 

\\ 






1 

1 

1 

1 


1 i 

i 

i 

/ 






w 

\ w 

, '- \ 
\ '■ \ 






1 

/ 

1 

/ 

1 

1 

/ 1 


y 

/ 

/ 

/ 






V • \ 


\ 

\ \ 

w 

• N \ 

. * . V 








/ 

/ 

/ 

/ 

/ 


1 i 
/ 

/ 

/ 

/ 












^ T~ — _ 






s 









-4 -2 2 4 6 8 10 12 

K 

FIG. 2: FCS generating function normalized to 2tTs//i at 
fixed temperature T/Tb = 0.5 and different voltages. 



merically for a wide range of parameters. In the second 
step one evaluates the P^s and calculates the cumulant 
generating function. 

Similarly to the free case the FCS turns out to be a uni- 
versal function of voltage V/Tb and temperature T/Tg 
measured in units of tunneling strength. With all funda- 
mental constants restored (spinful case) one obtains 

lnx(/0 = jT B F(eV/T B ,T/T B ) . (15) 

The function F(V/T B ,T/T B ) is plotted in Figs. ED and 
[2] for different temperatures and voltages. By virtue 
of normalization all curves go through the coordinate 
origin. However, all of them possess yet another zero 



point at precisely k = V/T = /3V, as is indeed pre- 
scribed by Eq. ([2]). Even more interesting is the overall 
curve symmetry with respect to reflections around the 
line k = V/2T, which is again perfectly consistent with 
the above fluctuation theorem. While this is not surpris- 
ing in the regimes of high temperatures and low voltages, 
when the system can be considered as being more 'clas- 
sical', it is most interesting to verify that this property 
persists at strong coupling as well. At T = the sec- 
ond zero moves to infinity and the curve approaches the 
'hockey stick' shape described by the analytical solution 
presented in [l2| , see Fig. [TJ 

To conclude, we presented a method to determine the 
full cumulant generating function of the nonequilibrium 
BSG model. It uses a modification of the thermodynamic 
Bethe ansatz and allows access to all values of the count- 
ing field and system parameters. At finite temperature 
we confirm the validity of the CGFT in all regimes and 
reproduce all known limiting cases. While the ultimate 
proof of the CGFT for generic interacting quantum sys- 
tem is still an open issue we have supplied an example 
of a nontrivial genuinely interacting system, in which it 
unmistakably holds. One avenue of future research might 
be the numerical evaluation of the FCS in the spirit of 
15|. 
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